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Abstract 

In this note we put forward a conjecture on the average optimal length for bipartite 
matching with a finite number of elements where the different lengths are independent 
one from the others and have an exponential distribution. 



The problem of random bipartite matching (or assignment) is interesting both from the 
point of view of optimisation theory and of statistical mechanics [1, 2, 3, 4]. 

Each instance of the problem is characterised by a N cross N matrix d; sometimes 
d(i, k) represents the distance between i and k. We are interested to compute the length 
of the optimal bipartite matching denned as: 

L(d)= mm d(i,U(i)), (1) 

i=l,N 

where the minimum is done over all the N\ permutations (II) of iV elements and IT(z) 
denote the results of the action of the permutation IT on z. 

If the matrix d has a given probability distribution, we can define the average of L{d) 
over this probability distribution. The problem of computing this average can be studied 
using the technique of statistical mechanics. 

A well studied case is when the elements of d are statistically independent one from 
the others. Here we can use replica theory [5, 6] in the mean field approximation. When 
the probability distribution of each element of d is flat in the interval [0, 1], after a long 
computation it was found [2, 3] that for large N we have: 

CP) + 2<(3) /J_\ 



(L)„ = C(2)- " ' N -' +O^J, (2) 

where ( is Riemann zeta function and the brackets denote the average over different in- 
stances of the system. This behaviour is well confirmed by numerical simulations [3, 4]. 
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In this note we consider a slightly different case where the probability distribution of 
each element of d is given by 

P(d) =exp(-d). (3) 

The result in the infinite N limit depends only on the behaviour of P(d) around d — 0, 
so that also in this case we expect that 

lim (L) N = C(2). (4) 

N—>oo 

Here we are interested in an exact formula for (L)n. Generally speaking (L)n is a 
rational number. A simple computation shows that 

(L) 1 = l (L) a = 1 + 1 (5) 

Our conjecture is that 

(L)n = E 4 (6) 

k=l,N ^ 

It is evident that the exact results for N — 1 and N = 2 and the limit for iV — > oo (eq. 
4) are compatible with our conjecture. 

We have tried to falsify numerically our conjecture for N = 3, 4, 5 (by randomly sorting 
O(10 8 ) instances of the matrix d). We have verified that the numerical results are com- 
patible with our conjecture with a precision of less that 5 10 -5 . It is remarkable that for 
these value of iV we have (L 2 (d)) ~ (L(d)) + 1. 

If this conjecture (eq. 6) is correct, it seems likely that the value (L) n can be computed 
analytically in a compact way. A proof of this result would be valuable also because it 
would automatically imply the correctness of the replica result (eq. 2) in the infinite N 
limit. However we have not succeeded in this task. 
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